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Abstract 

A new approach, which is based on the new canonical equations of Hamilton found by us recently, 
is presented to analytically obtain the approximate solution of the nonlocal nonlinear Schrodinger 
equation (NNLSE). The approximate analytical soliton solution of the NNLSE can be obtained, and 
the stability of the soliton can be analytically analysed in the simple way as well, all of which are 
consistent with the results published earlier. Eor the single light-envelope propagated in nonlocal 
nonlinear media modeled by the NNLSE, the Hamiltonian of the system can be constructed, which 
is the sum of the generalized kinetic energy and the generalized potential. The extreme point of the 
generalized potential corresponds to the soliton solution of the NNLSE. The soliton is stable when 
the generalized potential has the minimum, and unstable otherwise. In addition, the rigorous proof 
of the equivalency between the NNLSE and the Euler-Lagrange equation is given on the premise 
of the response function with even symmetry. 
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I. INTRODUCTION 


The propagations of the (l+D)-dimensional light-envelopes in nonlinear media have been 
stndied extensively for a few decades [l-7|, which are governed by the following dimensionless 
model js]: 

f)tn 

( 1 ) 


I^rup = 0 , 
oz 


where (p(r, z) is the complex amplitnde envelop, An(r, z) is the light-indnced nonlinear re¬ 
fractive index, 2 ; is the longitndinal coordinate, r is the U-dimensional transverse coordinate 
vector with D being the positive integer, and V_l is the U-dimensional differential operator 
vector of the transverse coordinates. Generally, A?7,(r, z) can be phenomenologically ex¬ 
pressed as the convolntion between the response fnnction R{r) of the media and the modnlns 
sqnare of the light-envelope (p{r, z) for the bnlk media with the nonlocal nonlinearity 8l-lll| 


A?7,(r, 2 ;)= / R{r — r')\(p{r\ z)\‘^d^r' 


( 2 ) 


According to the relative scale of the characteristic length of the response fnnction R and the 
scale in the transverse dimension occupied by the light-envelope (^, the degree of nonlocality 
can be divided into four categories 8|, llO, : local, weakly nonlocal, generally nonlocal, and 
strongly nonlocal, and locality is the case when the response function R is the Dirac delta 
function. In the local case, Eq.([I]) is reduced to 

.dip 
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( 3 ) 


Eq. (IT]) together with the nonlocal nonlinearity (|2|) is called as the nonlocal nonlinear 
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j while its special case, Eq. 

3 |. 


), is the well-known 


The NNLSE 
optical beams 


Schrodinger equation (NNLSE) 
nonlinear Schrodinger equation (NLSE) 

with its special case NLSE) can describe the nonlinear pro pag ations of the 
Oj, the optical pulses l|, yf and the optical pulsed beams |3|, l7|. The second 
term of the NNLSE accounts for the diffraction for the first case where r is the spatial 
transverse coordinate, the group velocity dispersion (GVD) for the second case where r is 
the time coordinate, and both the diffraction and the GVD for the last case where r is 
both the spatial transverse coordinate and the time coordinate, while the third term (the 
nonlinear term) describes the compression of the light-envelopes for all cases. Specifically, 


when D = 1, the NNLSE can model the propagation of the optical beam 


101 in the 
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self-focusing nonlinear planar waveguide, and can also model the propagation of the optical 


pulse [l| in the self-focusing nonlinear waveguide if the carrier frequency is in the anomalous 
GVD regime or in the self-defocusing nonlinear waveguide when its carrier frequency is in the 
normal GVD regime. The (1-1-1)-dimensional NNLSE has the spatial (or temporal) bright 
optical soliton solution When D = 2, the NNLSE can only describe the propagation 
of the optical beam in the nonlinear bulk media. The bright spatial optical soliton can 
exist stably for the nonlocal case [n|, but for the local case the strong self-focusing of a 
two dimensional beam will lead to the catastrophic phenomenon 12|. When D = 3, the 
NNLSE can describe the propagation of the optical pulsed beams. Like the case of D = 2, 
the self trapped optical pulsed beam propagating in the local nonlinear media will lead to 


the spatiotemporal collapse 


13| . which can be arrested by the nonlocal nonlinearity 


B- 


But when D > 3, the NNLSE is just a phenomenological model, the counterpart of which 
can not be found in physics. It’s important to note thatj^ the response function R is 
symmetric for the spatial nonlocality, but is asymmetric for the temporal nonlocality due to 
the causality jl^ . 

As the special case of the NNLSE, the NLSE ([3]) can be solved exactly using inverse¬ 
scattering technique j^, 16| when D = 1. But for the general case, a closed-form solution 
of NNLSE ([1]) cannot been found except for the strongly nonlocal limit, where the NNLSE 
can be simplihed to the (linear) Snyder-Mitchell model for the spatial nonlocality and an 
exact Gaussian-shaped stationary solution known as accessible soliton was found fl. Ap¬ 
proximately analytical solutions can be obtained by various of perturbation methods, such 
as the perturbation approach based on the inverse scattering transform 17|], the adiabatic 


perturbation approach 


variational method 


20 


18], the method of moments |^, and the most widely used one is 
23| . It was claimed without proof that the variationa 


method can 


only be applied in nonlocal cases where the response function is symmetric 2^. But the 
equivalency between the NNLSE ([T]) and the Euler-Lagrange equation is not proved rigor¬ 
ously until the mathematical proof given in the paper on the premise of the response function 
with even symmetry. And for the case of the response function without even symmetry, the 
method of moments can work well. Another new approach is presented here, and we apply 
the canonical equations of Hamilton to study the nonlinear light-envelope propagations. By 
taking this approach, the approximate analytical soliton solution of the NNLSE is obtained. 
Furthermore, the stability of solutions can be analysed analytically in a simple way as well. 
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but it can not be done by the variational approach. 

The paper is organized as follows. We hrstly give the rigorous proof of the equivalency 
between the NNLSE and the Euler-Lagrange equation in Sec. [ITl which is the basis of the 
variational approach applied in the NNLSE. The canonical equations of Hamilton (CEH) 
is a parallel method to the Euler-Lagrange equation in classical mechanics. But we hnd 
that the conventional CEH can not restate the NNLSE, and present a new CEH to restate 
the NNLSE, which is outlined in Sec. m Based on the new CEH, we introduce a new 
approach in Sec. [IV] to deal with the nonlinear light-envelope propagations. In Sec. |V] two 
remarks on the new approach are made. Firstly, we show that the conventional CEH will 
yield contradictory and inconsistent results. Secondly, we discuss the differences between 
our approach and the variational approach. Sec. IVII gives the summary. 


II. PROOF OF THE EQUIVALENCY BETWEEN THE NONLOCAL NONLIN¬ 
EAR SCHRODINGER EQUATION AND THE EULER-LAGRANGE EQUATION 


The variational approach is a widely used method to obtain the approximately analytical 
solution of the NLSE 
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21| . The reason why the variational approach can be used is that 


the NLSE can be restated by the Euler-Lagrange equation, which reads (for the sake of 
simpleness, only the case that D = 1 is taken consideration here) 


d dl d dl 
Wzd{%-) 


if the Lagrangian density I is given by 201 
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(4) 


(5) 


Replacing ip* with ip, the complex-conjugate equation of the NLSE can be obtained from 
the Euler-Lagrange equation (jl]) consistently. Although the variational approach has been 
applied to the problems associated with the NNLSE ([1]), in which the Lagrangian density is 
expressed as |22|, [2^ 


with Ini = ||(p(r, 2 :)pAn(r, z), the equivalency between the NNLSE ([1]) and the Euler- 
Lagrange equation has not been proved rigorously. Without loss of generality, here we only 
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give the proof of the equivalency in the case that D = 1, and the extension to the general 
case of D will be easy in a similar way. 

Comparing the two expressions of the Lagrangian density for the NLSE and the NNLSE, 
i.e., Eqs. (E]) and ([H]), we can observe that the Lagrangian density for the NNLSE contains a 
convolution between the response function R and the modulus square of the light-envelope 
(f. Therefore, it has been somewhat confused how to calculate such terms as dlni/dip* and 
dlni/d(p for the NNLSE since Ini is not the function of (p and ip* but the functional of them. 
To this end, we hrst construct a functional as 

/ oo poo 

/ Inidxdz 

■OO J —OO 


1 
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f^oo roo roo 


R{x — x^)\lp{x\ z)\ \^{x^ z)\^dx'dxdz. 


( 7 ) 


— OO ^ —OO ^ —oo 


The variation of the functional F{p, p*) is dehned as 25| 

SF{p,p*) = -^Fi^p + e6p,p* + e6p*)\e^o 


1 
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^oo roo roo 


R{x — x')\p{x^ z)\^ z)6p*{x\ z) -\- p*{x\ z)6p{x\ z)] dx'dxdz 


— OO J —OO J —OO 

roo roo 


' —oo J —oo 


An [p{x, z)5p*{x, z) p*{x, z)5p{x, z)] dxdz. 

If the response function is symmetric, i.e., R{x) = R{—x), then we can obtain that 
R{x — x')\p{x, z)\‘^ i z)dp*{x', z) + p*{x', z)5p{x', z)] dx'dx 

= I An[p{x, z)6p*{x, z) + p*{x, z)6p{x, z)] dx. 

J —OO 

Then the variation of the functional F{p, p*) is simplihed to 


( 8 ) 
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( 9 ) 


poo roo 


5F{p,p*) = / Anp*{x^z)5p{x^z)dxdz R 


Anp{x, z)6p* {x, z)dxdz. (10) 


' —OO J —oo 


On the other hand, the variation of the functional F{p, p*) can be also expressed as 26| 


POO roo 


6F{p,p*) = 6 


lni(,P^,P^*)dxdz 


’ —oo J — OO 
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dlni 

d dlnl 

d 

dlnl 

J —oo J —oo 

dp 


dx 
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1_ 

d dlnl 

1 

1- 
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dp* 


dx 
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Spdxdz 

6p*dxdz 
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1-00 dp 


Spdxdz + 
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nl 


1-00 dp^ 


5p*dxdz. 


( 11 ) 


5 
















Comparing Eqs. ffT0|) and ffTTj) . we obtain 


dl 


nl 


dcp* 

dlnl 

dip 


= Anip{x,z), 
= An(p*{x, z). 


( 12 ) 

(13) 


Inserting the Lagrangian density (E]) for the case of D = 1 into the Euler-Lagrange equation 
(jH), the hrst two terms of (jl]) can be easily obtained as 

d dl d dl 


.dip ^ d‘^ip 


(14) 


dx d d dz dx"^' 

Then the NNLSE ([1]) can be obtained from the Euler-Lagrange equation (|1]) by using 
Eqs. m and ffTTD respectively, and its complex-conjugate equation can also be obtained 
consistently in a similar way. 

Consequently, the NNLSE ([T]) is equivalent to the Euler-Lagrange equation (jl]) if the 
response function is symmetric. But for the asymmetric response function, for example, the 


response function for the temporal nonlocality 1^, we can not show the equivalency between 
the NNLSE ([1]) and the Euler-Lagrange equation (jlD anymore. In our points of view, the 


conclusion obtained here is equivalent to that given in Ref. 2^, where the authors claimed 
that the equation (1) in Ref. 2^ (similar to the NNLSE) does not have a Lagrangian when 
the temporally asymmetric nonlocal term is included and that “Had the nonlocality been 
symmetric, then variational techniques could have been applied”, although no any proof was 
given in Ref. 241. 


III. CANONICAL EQUATIONS OF HAMILTON FOR THE NNLSE 

As discussed in the section above, the variational approach to hnd the approximately 
analytical solution of the NNLSE is based on the Euler-Lagrange equations. In the classical 
mechanics, however, there exist two theory frameworks: the Lagrangian formulation (the 
Euler-Lagrange equations) and the Hamiltonian formulation (canonical equations of Hamil¬ 
ton). The two methods are parallel, and no one is particularly superior to the another for 


the direct solution of mechanical problems 


27l |. The new approach presented in this paper 


to analytically obtain the approximate solution of the NNLSE is based on the new canonical 
equations of Hamilton (CEH) found by us recently 28|. For the sake of the systematicness 
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and the readability of this paper, the key points abont the new CEH are ontlined here in 


this section, althongh the detail can be fonnd in Ref 


28|. 


We hrstly dehne two different systems of mathematical physics: the second-order differ¬ 
ential system (SODS) and the hrst-order differential system (FOBS). The SODS is dehned 
as the system described by the second-order partial differential eqnation abont the evolntion 
coordinate, while the FOBS is dehned as the system described by the hrst-order partial dif¬ 
ferential eqnation abont the evolntion coordinate. The Newton’s second law of motion and 
the NNLSE are the exemplary SODS and FOBS, where the evolntion coordinates are the 
time coordinate t and the propagation coordinate 2 ;, respectively. The conventional CEH 

dH 

dpi' 

dH 
dqi ’ 

is established on the basis of the Newton’s second law of motion. The dot above the variable 
in Eqs. flT^ and flTbl) (g* and Pi) indicates the derivative with respect to the evolntion coor¬ 
dinate (here the evolntion coordinate is the time t), Qi and pi are said to be the generalized 
coordinate and the generalized momentnm, and H is the Hamiltonian. The CEH (|T5|l and 
ffTbjl can be extended to the continnons system as |27| 


Qi = 

-Pi = 


(15) 

(16) 


Qs = 


6h 

SPs' 

6h 

-Ps = 


( 17 ) 

(18) 


with s = 1, • • • , iV representing the components of the qnantity of the continnons system |2 


Sh 

Sqs 


dh 

dqs 


and 4^ = denote the fnnctional derivatives of h with respect 

to qs and ps with qs^x = ^ and Ps^x = and h is the Hamiltonian density of the continnons 

system. 

We have shown that the FOBS can not be expressed by the conventional CEH, and we 
have re-constrncted a set of new CEH throngh the following procednre. 

For the hrst-order differential system of the continnons systems, the Lagrangian density 
mnst be the linear fnnction of the generalized velocities, and expressed as 


N 


^ ^ ^ Rs(gs)gs T Q(^qsy qs,x)^ 


(19) 


S=1 
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where Rs is not the function of a set of with = dqs/dx. Consequently, the generalized 
momentum which is obtained by the dehnition ps = dl/dqs as 


P, = = ,N), 


( 20 ) 


is only a function of There are 2N variables, g^ and ps, in Eqs. fl20|) . The number of Eqs. 
(|20|) is N, which also means there exist N constraints between g^ and Ps- So the degree of 
freedom of the system given by Eqs. (1^ is N. Without loss of generality, we take gi, • • • ,qu 
and Pi, • • • iPfi as the independent variables, where v + p = N. The remaining generalized 
coordinates and generalized momenta can be expressed with these independent variables as 
qa = ga(gi, • • • ,g, 2 ,Pi, • • • ,p^x)(a = I/ + l, • • • ,iV), and P/? = pp{qi, ■ ■ ■ ,g^,pi, • • • ,p^,){|3 = 
p + 1, • • • iN)- The Hamiltonian density h for the continuous system is obtained by the 
Legendre transformation as h = qsPs — J-, where the Hamiltonian density h is a function 
of u generalized coordinates, gi,--- ,qu, and p generalized momenta, pi,-- - ,p^. We can 
obtain the new CEH consisting of N equations as 



( 21 ) 



( 22 ) 


where A = 1, • • • , i/, p = 1, • • • , p, and v + p = N. The CEH (1^ and (12^ can be easily 
extended to the discrete system, which can be expressed as 



(23) 



(24) 


where A = 1, • • • , i/, p = 1, • • • , p, v R p = N^ the generalized momenta are dehned as 


dL 

Pi = 


(25) 


with L = Id^r being the Lagrangian, and the Hamiltonian is obtained by Legendre 
transformation as 


n 
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( 26 ) 











For the SODS, all the generalized coordinates and the generalized momenta are independent, 
the new CEH fl2^ and fl2^ are automatically reduced to the conventional CEH flT^ and 


(uni). 


We have shown that the FOBS can only be expressed by the new CEH, but do not by 
the conventional CEH, while the SODS can be done by both the new and the conventional 
CEHs. We have also shown that the NLSE can be expressed by the new CEH in a consistent 
way if the propagation coordinate ^ in the NLSE is considered to be the evolution coordinate. 

IV. APPLICATION OF THE NEW CEH TO LIGHT-ENVELOPE PROPAGA¬ 
TIONS 

Different from the case of the NLSE, the Hamiltonian density of the NNLSE contains 
the convolution between the response function R{r) and the modulus square of the light- 
envelope ip{r, z). Following the procedure in Sec. [TTl h can be easily proved that the NNLSE 
can also be expressed with the new CEH in a consistent way if the propagation coordinate z 
in the model is considered to be the evolution coordinate. Based on the new CEH, we now 
introduce a new approach to deal with the nonlinear light-envelope propagations. 

We assume the trial solution of the form as 



( 27 ) 


where qa, qe are the amplitude and phase of the complex amplitude of the light-envelope, 
respectively, is the width of the light-envelope, qc is the phase-front curvature, and they 
all vary with the propagation distance (the evolution coordinate) ; 2 . The response function 
of materials is assumed as 



( 28 ) 


Inserting the trial solution (HT} and the response function fl28l) into the Lagrangian density 
(1^ . and performing the integration L = we obtain 



( 29 ) 


which is a function of generalized coordinates, qA, qw, qc and generalized velocities, gc, qe (The 


dot above the variable indicates the derivative with respect to the evolution coordinate z), 


9 






but not an explicit function of the evolution coordinate z. Eq. fl^ can be understood as 
a “coordinate transformation”. Through such a transformation (of course, this is not a 
real coordinate transformation in the rigorous sense in mathematics), the coordinate system 
consist of a set of generalized coordinates (p and p* is transformed to that consist of another 
set of generalized coordinates gc, and g^, and the Lagrangian density expressed by 

Eq. (j6]) in the continuous system is transferred to the Lagrangian expressed by Eq. (1291) in 
the discrete system at the same time via the integration L = Id^r. 

Then the generalized momenta can be obtained by dehnition fl25D as follows 


Pa = Pw = 0, (30) 

p, = (31) 

P» = - (^) ‘ IaiP (32) 

The Hamiltonian of the system then can be determined by Legendre transformation fl2B]) 

+2^"^ 2-0(^7^ + g^)'^/^(l + g^gt)]) (33) 


and can be proved to be a constant, i.e. H = 0. 

There are four generalized coordinates and four generalized momenta in the four equations 
(130]) . fl3T]) and (133]) . So the degree of freedom of the set of equations (130]) . (I3T]) and (133]) 
is four. Without loss of generality, we take qc,<le,Pc and pe as the independent variables. 
By solving Eqs. (13T]) and (]33]) . the generalized coordinates qa and g^, can be expressed by 
generalized momenta Pc and pe as qa = {—peY^'^[Dp 0 /{ 2 'Kpc)]^^'^ and g^ = [4pc/(T*Pe)]^^^, 
and inserting this result into the Hamiltonian (133]) yields 

^ _D‘^pI + IQplql _ 1 D/2, 2 yD/2 ^343 

By use of the canonical equations of Hamilton (1331) and (133]) in the case that p = v = 2 and 
n = A because there are only two independent generalized coordinates and two independent 
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( 35 ) 


generalized momenta, we can obtain the following four equations 

, D-‘pI , , , 

^pI ^Pc + Dpewl^ 

. _ _{4^ +D)7r-^/‘^p^pe{^+ 

^Pc + Dpewl^ 

P,2 n -D/2 2 2 ( ipc I 2 \-D/2 

2pc 4pc + Dpswl^ 

Pc = 8pc()c, 

P9 = 0. 


(36) 

(37) 

(38) 


It can be found that the generalized coordinate qe is not contained in the Hamiltonian 
then qq is a cyclic coordinate. It is known that the generalized momentum conjugate 
to a cyclic coordinate is conserved j^. Therefore, the generalized momentum pg conjugate 
to the generalized coordinate qe is a constant, which can be conhrmed by Eq. (155]) . In fact, 
this represents that the power of the light-envelope, 

/ OO 

W? = q\{\/^qu,)^, (39) 

-OO 

is conservative. Then we can obtain 


9^ = a(vT2«») °' 


(40) 


Taking the derivative with respect to z on both sides of Eq.(|3T]), then comparing it with 
Eq.(|37]), we can obtain with the aid of Eq. fHOj) 


Qc = 




Then by substituting Eq. fHip into the Hamiltonian 
H = T + V, where 


(41) 

with the aid of Eq. (14011 . we have 


T = 4dP„4/, (42) 

V = K + -‘>1^ (43) 

qw ^ 

are the generalized kinetic energy and the generalized potential of the Hamiltonian system, 
respectively. 

Now we can observe that the dynamics of the light-envelopes in nonlinear media can be 
treated as problems of small oscillations of a Hamiltonian system about positions of equilib¬ 
rium from the Hamiltonian point of view. The equilibrium state of the system described by 
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the Hamiltonian given together by Eqs. fH2|) and (143|) corresponds to the soliton solutions of 
the NNLSE, and can be obtained as the extremum points of the generalized potential of the 
Hamiltonian system. An equilibrium position is classihed as stable if a small disturbance 
of the system from equilibrium results in small bounded motion about the rest position. 
The equilibrium is unstable if an inhnitesimal disturbance eventually produces unbounded 
motion 27|]. It can be readily seen that when the extremum of the generalized potential is 


a minimum the equilibrium must be stable, otherwise, the equilibrium must be unstable. 
In this sense, therefore, the viewpoint in some literatures 29l-l32| . where solitons were re¬ 


garded as the extremum of the Hamiltonian itself rather than the generalized potential 
the Hamiltonian system, would be some ambiguous. Because in those literatures 29 - 


of 

33 


the trial solution has a changeless prohle (solitonic prohle), the state expressed with the 
solitonic prohle is the static state. The kinetic energy of the static state is zero, and the 
Hamiltonian is equal to the potential of the static state. In this connection, the extremum 


of the Hamiltonian equals to the extremum of the generalized potential of t 


only in value. Although the soliton solutions obtained in such literatures 29 


le static system 


32| are correct. 


it is more reasonable to consider the soliton solutions of the NNLSE as the extremum points 
of the generalized potential of the Hamiltonian system. 

In order to hnd the equilibrium position (the soliton solution), letting dV/dqw = 0, we 
have 


32 


- — + Stt ^^^PoQu, {wl, + ql) 




= 0 . 




(44) 


We can easily obtain the critical power 

4 ^°/=« + 


Po = 


(45) 


with which the light-envelope will propagate with a changeless shape. In addition, when 
Pq = Pc, it can be easily obtained that (jc = qc = 0, which implies that the wavefront of the 
soliton solution is a plane. 

Then we elucidate the stability characteristics of the soliton by analysing the properties 
of the generalized potential V. Performing the second-order derivative of the generalized 
potential V with respect to q^,, then inserting the critical power into it, we obtain 


d^V 


64 

r. 2+p 1 

dql 

II 

1 

•45 1 

r 2(l + a2)J 


( 46 ) 
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where a = Wm/Qw is the degree of nonlocality. The larger is a, the stronger is the degree 
of nonlocality. When T > 0, the generalized potential has a minimum, and the soliton is 
stable. From Eq. fl46p we can obtain the criterion for the stability of solitons, that is 




(47) 


which is, in fact, consistent with the Vakhitov-Kolokolov (VK) criterion 33|(for detail, see 
the footnote 341). 


The local case 


When Wm —)■ 0, the response function R{r) —)■ 5(r), then the NNLSE will be reduced to 
the NLSE ([3]). In this case, Eqs. fjTHl) and fH6l) are reduced to 


32 


P, = 4;^°''y“-^T = -(2-C). 


(48) 


When D = 1, the critical power is deduced to Pc = which is consistent with Eq.(42) 


of Ref. 


20l |. When D = 2, the critical power is deduced to Pc = dvr, which is the same as 


Eq.(16a) of Ref. 36|. We can obtain T > 0 when D < 2, T < 0 when D > 2, and T = 0 
when D = 2. So for the local case, the soliton is stable for (l+l)-dimensional case, but 
unstable when D > 2. It needs the further analysis for the case of D = 2 because T = 0. 
When D = 2, the generalized potential fl4^ from the Hamiltonian point of view is deduced 
to 


R = 


(dTT - Pq) Pq 


(49) 


which has no extreme when Pq ^ dvr. When Pq = Pc = dvr, it can be obtained that V = 0, 
which is the extreme but not the minimum. So the (l+2)-dimensional local solitons are 
unstable. The relation between the potential V and the width qcu of the light-envelope is 
shown in Fig{T] If the power of the light-envelope equals to the critical power, the potential 
will be a constant, as can be seen by dash curve of Fig. ([1]). Without the external disturbance, 
the light-envelope will stay in its initial state, and keep its width changeless. But the ideal 
condition without external disturbances can not exist in fact. If the external disturbance 
makes the power larger than the critical power, then the system will evolve towards the 
lower potential, the beam width will become more and more smaller, and the optical beam 
will collapse at last, as can be conhrmed by the dash-dot curve of Fig{Tl If the external 
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disturbance makes the power smaller than the critical power, then the system will also 
evolve towards the lower potential, the beam width will become more and more larger, and 
the optical beam will diffract at last, as can be conhrmed by the solid curve of FigJT] These 


conclusions are consist with those of Refs. 


37 


39|. 



FIG. 1: Qualitative plot of the potential F as a function of qw for three cases, Pq < Pc (solid 
curve), Pq = Pc (dashed curve), and Pq > Pc (dash-dot curve) when D = 2. 


B. The nonlocal case 


For the nonlocal case, when D < 2, the condition flTTl) can be satished automatically. 
That is to say the (l+l)-dimensional and the (l+2)-dimensional nonlocal solitons are always 
stable when the response function of the material is a Gaussian function. It is consistent 
with the conclusion of Ref. 3^. When D > 2 the solitons can be stable only if the the 


degree of nonlocality is strong enough that can satisfy the criterion flTTl) . which is also the 
same as the result of Ref. [35 1. 


V. TWO REMARKS 

At the end, we make two remarks on the new approach in dealing with the nonlinear 
light-envelope propagations presented in the paper. Firstly, the new approach is based on 
the new CEH fl23|) and fl2T)) . and we will show that the conventional CEH (IT^ and flT^ will 
yield contradictory and inconsistent results. Secondly, we will compare our approach with 
the variational approach, and discuss the differences between them. 
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A. Contradictory results coming from the conventional CEH 


Here we use the conventional CEH ffT^ and ffT6|) to deal with the light-envelope prop¬ 
agated in nonlinear media, following the same procedure in Sec. IIVI and show that the 
conventional CEH 0151) and 016p will give the contradictory and inconsistent results. 

Without loss of generality, we only take the NLSE (|3]) with H = 1 as an example. The 
NLSE is a special case of the NNLSE when Wm approaches to zero. Then letting Wm = 0 
and D = 1 makes Hamiltonian fl33|) reduced into 


[2V2 - ql {q\ - 2 ^qlqfi] 

" =- 

Because the Hamiltonian is only the function of the generalized coordinates, the CEH ffl^ . 
the right hand side of which is the derivative of the Hamiltonian with respect to the gen¬ 
eralized momentum, can yield nothing unless (jc = Qe = Qa = Qw = 0. It means the four 
quantities are all the conserved quantities. This result coming from the CEH flT5|l is obvi¬ 
ously wrong because such quantities as the amplitude qA, the width q^, and the phase-front 
curvature qc all generally vary with the evolution coordinate z except for the soliton state, 
and qg, the phase of the complex amplitude of the light-envelope, must be the function of 2 ; 
even for the soliton state. 

From the other CEH ffT6|) . four equations can be obtained as 


Pc 

Pe 

Pa 

Pw 


dH 

dqc 

dH 

dqe 

dH 

dqA 

dH 

dq-u, 


-V^qiqcql, 


0 , 

y/^qA {-V 2 + q\ql - \/2ggg^) 

Qw 

- sAight) 
Mi 


(51) 

(52) 

(53) 

(54) 


Substitution of the generalized momenta Pc given by Eq. fl3ip into Eq. (IHip yields the 
same result as Eq. (HI]). Then inserting Eq. flTip into the Hamiltonian (l50P gives out 
H = ^ (^q'^ + ^ — ^ . The Hamiltonian is the sum of the generalized kinetic energy 

and the generalized potential V{qw) = ^ (4- 'p~): which is also the same as Eq. (Hsp 

when D = 1 and Wm = 0. Therefore, the critical power, corresponding to the extremum 
point of the generalized potential. Pc = is the same as Eq. fHsp when D = 1 and Wm = 0. 
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It can also be found that Eq. fl5^ is the same as Eq. fl38D . which means that the power of 
the light-envelope is conservative. Although the hrst two equations, Eqs. flSTj) and ([52]), of 
a set of equations resulting from CEH flT^ can give out the correct results, the other two 
equations, Eqs. and fl51|l . will yield the contradictory and inconsistent results. Let us 
show as follows. Inserting Eq. (150]) into Eqs. (105]) and (105]) yields 

S-v/vr 


Pn = 


Qc = 


5qw ’ 

3 1 

5 ^ 


(55) 

(56) 


Obviously, the two results given by Eqs. (100]) and (106]) are both wrong. Under the assumption 
of the light-envelope with the form of Gaussian-shape given by Eq. fl27]) . the power carried 
by the light-envelope should be Pq = \/t^/‘ 2^q\qw given by Eq. f]59]) . with which Eq. f]55]) is 
contradictory and inconsistent. Eq. flOO]) gives the hxed relation between and But the 
phase-front curvature, gc, should be changed depending upon the state of the light-envelope, 
especially qc should be zero for the soliton state, with which Eq. (106]) is inconsistent. 

It is no surprise to obtain such contradictory and inconsistent results from the canonical 
equations of Hamilton (ITO]) and (1T6]) . since both the NNLSE ([1]) and its complex conjugation 
can not be derived from the canonical equations of Hamilton (1T7]) and (1T8]) as stated in 

Sec. [ml 


B. Our approach vs the variational approach: same and different 


As mentioned above, our approach presented in this paper is based on the canonical 
equations of Hamilton (the Hamiltonian formulation), while the variational approach 20|] is 
based on the Euler-Lagrange equations (the Lagrangian formulation). Although the same 
point of the two approaches is to hrst compute the Lagrangian of the system by using a 
suitably chosen trial function, they are in essence two parallel methods because the Hamil¬ 
tonian formulation and the Lagrangian formulation are two parallel theory frameworks in 
the classical mechanics. 

The most important concept in our approach is the “potential”. The potential given 
by Eq. f]55]) is the real “potential” of the system that a single light-envelope propagates in 
nonlocal nonlinear media modeled by the NNLSE. It is not, of course, the potential of the 
narrow-sense mechanical system, but does be the potential in the frame of the Hamiltonian 
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theory, that is, the potential of the Hamiltonian system. In other word, it is the potential 
from the Hamiltonian point of view. Looking back to the variational approach, we can ob¬ 
serve that although the “potential” was also introduced [see, Eqs. (28) and (29) in Ref. j^]. 
it is just a mathematically equivalent potential in the sence that the evolution of the width 
of the light-envelope can be analogous to that of a particle in a potential well, rather than 
the real “potential” of the system. 


VI. CONCLUSION 

We introduce a new approach, based on the new canonical equations of Hamilton found 
by us recently, to analytically obtain the approximate solution of the nonlocal nonlinear 
Schrodinger equation and to analytically discuss the stability of the soliton. For the single 
light-envelope propagated in nonlocal nonlinear media modeled by the NNLSE, the Hamil¬ 
tonian of the system can be constructed as the sum of the generalized kinetic energy and 
the generalized potential. The extreme point of the generalized potential corresponds to the 
soliton solution of the NNLSE. The soliton is stable when the generalized potential has the 
minimum, and unstable otherwise. In addition, we give the rigorous proof of the equiva¬ 
lency between the NNLSE and the Euler-Lagrange equation on the premise of the response 
function with even symmetry. 
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